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Abstract 

Events with a (2+1) jet topology in deep-inelastic scattering at HERA are studied in the 
kinematic range 200 < < 10000 GeV^. The rate of (2+1) jet events has been deter- 
mined with the modified JADE jet algorithm as a function of the jet resolution parameter 
and is compared with the predictions of Monte Carlo models. In addition, the event rate is 
corrected for both hadronization and detector effects and is compared with next-to-leading 
order QCD calculations. A value of the strong coupling constant of q;s(M^)=0.118 it 
0.002 (stat.) Iqoos (■sys*-) ooe {theory) is extracted. The systematic error includes un- 
certainties in the calorimeter energy calibration, in the description of the data by current 
Monte Carlo models, and in the knowledge of the parton densities. The theoretical error is 
dominated by the renormalization scale ambiguity. 
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1 Introduction 



We present a study of events with a (2+1) jet topology and determine the strong coupling con- 
stant, as, using neutral current deep-inelastic scattering (DIS) events recorded with the HI 
detector at the ep collider HERA in 1994 and 1995. In this period HERA was operated with 
positron and proton beams of 27.5 and 820 GeV energy, respectively, corresponding to a centre- 
of-mass energy of ^/s = 300 GeV. 

In the Quark-Parton-Model, neutral current DIS corresponds to the interaction of a virtual 
photon or Z° boson with a quark in the proton. The interaction can be characterized by the 
two independent variables and x where is the absolute value of the virtual boson 4- 
momentum squared and x is related to the fraction of the proton momentum carried by the 
struck quark. Experimentally, events with a (1+1) jet topology are observed. The notation '+1' 
refers to the proton remnant jet. QCD corrections in 0{as), namely QCD-Compton scattering 
ilQ ^ Q9) ^iid Boson-Gluon-Fusion (7(7 qq), lead to (2+1) parton final states. Due to the 
high centre-of-mass energy at HERA, multi-jet structures have been observed clearly [jl]], and 
quantitative tests of QCD and the determination of the strong coupling constant are made 
possible. 

Previous jet analyses and determinations of at HERA were based on the measurement of 
R2+1 iQ'^), the (2+1) jet event rate as a function of [0. The jets were found by applying the 
modified JADE jet algorithm [j^] in the laboratory frame for a fixed value of the jet resolution 
parameter. In particular, the measurement of -R2+i(Q^) allows the dependence of on the 
scale to be studied in a single experiment. In this analysis a complementary approach is 
adopted. We take events in the range 200 < < 10 000 GeV^ and then cluster measured 
calorimeter energy depositions with the modified JADE jet algorithm until (2+1) jets remain. 
The minimum mass squared of any pair of the (2+1) jet four- vectors, scaled by the hadronic 
energy squared W'^, is the variable y2, which we study. For a clear (1+1) jet event a small value 
of y2 is expected whereas any event with a larger jet multiplicity must result in a large value. 

This is the first measurement of differential jet event rates at HERA [|p. The presence of 
a strongly interacting particle in the initial state gives rise to considerable differences from the 
situation in e+e annihilation, where differential jet event rates have been studied in much detail 
[Q]. The proton remnant, the initial state QCD radiation, the large momenta of the produced jets 
in the direction of the incoming proton, and finally the uncertainties in the knowledge of the 
parton content of the proton complicate the measurement. The study of the same observable in 
processes as different as e+e^ annihilation and ep scattering, however, may lead to improved 
understanding of systematic uncertainties in the determination of from hadronic final states 
and provides an important test of QCD. 

The analysis consists of the following steps. After the data selection, the accuracy with 
which the data are described by the colour dipole Monte Carlo model ARIADNE [§] and the 
leading-logarithm parton shower model LEPTO [[7]] is studied. The rate of (2+1) jet events 
is corrected for detector acceptance, resolution and inefficiencies as well as for hadronization 
effects. A sophisticated correction procedure is used that takes migration effects into account. 
Next, the parton jet distributions of these models are compared qualitatively with next-to- 
leading order (NLO) calculations available in the form of the programs MEPJET ^ and DIS- 
ENT , in order to verify that a jet phase space region has been selected in which the NLO 
calculations can be expected to be a good approximation to the data. Finally, the NLO calcu- 
lations are fitted to the corrected data as a function of a^, and the systematic uncertainties are 
evaluated. 
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2 The HI detector 



A detailed description of the HI detector can be found in [|T0|]. The components most relevant 
for this analysis are the central tracking system, the liquid argon calorimeter, the backward 
electromagnetic calorimeter, and the instrumented iron return yoke. 

The central tracking system consists of several inner and outer drift and proportional cham- 
bers. It is used in this analysis to determine the ep collision point and to aid the identification 
of the scattered positron. The tracking system is surrounded by a large liquid argon sampling 
calorimeter covering a polar angle range of 4° < 6' < 154°. The polar angle 6 is measured with 
respect to the incoming proton beam which is defined to point in the -^z direction. The electro- 
magnetic and hadronic sections of the liquid argon calorimeter correspond in total to a depth of 
4.5 to 8 interaction lengths. The energy resolution of the liquid argon calorimeter for electrons 
and hadronic showers is a/E = 12%/^E{GeV) © 1% and a/E = 50%/^E{GeV) © 2%, 
respectively [pT]]. The absolute energy scale for hadronic energy depositions is known to better 
than 4%, and that for electromagnetic energy depositions to better than 3%. 

Since 1995 the backward region of the HI detector has been equipped with a drift cham- 
ber and a lead/scintillating-fibre calorimeter. Its main purpose is the detection of electrons at 
small scattering angles. In addition, the timing information it provides allows efficient dis- 
crimination against out-of-time proton beam related background events at early trigger levels. 
Before 1995 the backward region was instrumented with a multi-wire proportional chamber, a 
lead/scintillator electromagnetic calorimeter, and a scintillator array for timing measurements. 

Outside the calorimeters a large superconducting solenoid provides a magnetic field of 
1.15 Tesla. The instrumented iron return yoke identifies energetic muons and detects leakage of 
hadronic showers. 



3 Event and cluster selection 

Neutral current DIS events are selected using the following criteria. We require a scattered 
positron candidate to be detected within 9e < 150° so that it is well contained within the ac- 
ceptance of the liquid argon calorimeter. A cluster of contiguous energy depositions in the 
calorimeter is identified as a positron candidate if its energy deposition in the electromagnetic 
calorimeter section exceeds 80% of the cluster energy and if its lateral and longitudinal profiles 
are compatible with those of an electromagnetic shower [[T2|]. In addition, its position must be 
matched to a reconstructed track to better than 1 .7° in polar angle and to better than 6° in az- 
imuthal angle. The range is restricted to 200 < < 10000 GeV^ where Q"^ is determined 
from the scattered positron energy and polar angle. The cut Q'^ > 200 GeV^ offers several 
advantages: hadronic final state particles are better contained in the detector since they must 
balance the transverse momentum of the scattered positron, which is detected in the liquid ar- 
gon calorimeter at large Q^; the range of x is implicitly restricted to larger values of x where 
the parton density of the proton is better known and where initial-state QCD radiation beyond 
NLO is suppressed. 

The measured z coordinate of the primary event vertex is required to be within a distance 
of 30 cm from the nominal ep collision point. The time-of-flight information of the backward 
scintillator array is required not to be inconsistent with impact times of particles originating 
from the ep collision point. Both cuts strongly reduce proton beam-related background events. 



5 



The inelasticity y = Q'^/sx, calculated from the scattered positron energy and polar angle, 
is required to be smaller than 0.7. This cut corresponds to a polar angle dependent minimum 
positron energy requirement to suppress background from misidentified photoproduction events 
and to reduce the influence of QED radiation. The remaining effects of initial and final state 
QED radiation were studied with DJANGO [[T3|]. They were found to be small and are neglected 
in the following. Photoproduction and beam-related background events are further suppressed 
by requiring 30 < E-P^ < 70 GeV where E and P^, are the summed energy and longitudinal 
momentum components of all reconstructed clusters, each assumed to be massless. For NC DIS 
events E-P^ is ideally expected to be 55 GeV, corresponding to twice the positron beam energy. 
The invariant mass squared of the hadronic final state, VF^i, as calculated using the double angle 
method [ [T4l ] is required to exceed 5 000 GeV^ to ensure a substantial hadronic activity for jet 
production. In addition, we reject events where cosmic muons or beam halo muons crossing the 
detector are identified [[T3|]. 

The events recorded were triggered by the electron trigger of the liquid argon calorimeter. 
The above cuts imply that the energy of the scattered positron always exceeds 10 GeV. The 
average trigger efficiency for the selected data sample was found to be larger than 99% and is 
independent of the hadronic final state. 

With these cuts we obtain a sample of 11 192 deep-inelastic scattering events correspond- 
ing to an integrated luminosity of 7 pb~^. The remaining background from beam-gas collision, 
photoproduction, cosmic muon or halo muon events in this sample is negligible. In particular, 
the fraction of photoproduction events is estimated to be less than 0.5%. The fraction of diffrac- 
tive events [|1^, defined as events with a calorimeter energy deposition of less than 0.5 GeV in 
a cone of 15° around the beam direction, is of the order of 1%. After the application of the jet 
algorithm, further cuts are applied to select a subsample enriched with (2+1) jet events. 

In this analysis, hadronic jets are reconstructed from the energy depositions in the liquid 
argon calorimeter and the instrumented iron. Clusters that are not well measured or that are 
not related to the hadronic final state are rejected by the following quality cuts: the polar an- 
gle, 9dus, of a cluster is required to satisfy Odus > 7° to select clusters that are well within the 
geometrical acceptance of the liquid argon calorimeter, and energy depositions in the backward 
electromagnetic calorimeter are discarded since this has limited containment for hadrons. Fur- 
ther requirements of less importance are: the energy fraction leaking into the instrumented iron 
is required to not exceed 40%; hadronic clusters must be separated from the positron candidate 
by an angle greater than 10°; clusters with an angle of larger than 50° with respect to their 
closest neighbouring cluster are rejected. This latter cut is imposed to decrease the sensitivity 
to isolated noise contributions or to photons radiated from the scattered positron. After these 
selections the average number of accepted clusters per event is 37.8. 

4 Jet algorithms and jet event rate definition 

The jets in a given event are found using the JADE jet algorithm [^. The jet algorithm is applied 
in the laboratory frame to the clusters of the liquid argon calorimeter and the instrumented iron 
satisfying the cuts given in section 3. The algorithm is modified compared to the version used 
in e+e~-annihilation in two respects: (a) the cluster that is attributed to the scattered positron is 
removed; (b) a massless four-vector is determined and is treated as an additional cluster by the 
jet algorithm to account for the longitudinal component of the momentum carried by the proton 
remnant particles escaping through the beam pipe. 
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The jet algorithm calculates the scaled quantity mjj/W'^ of pairs of clusters or 'proto' jets 
i, j, where W'^ is the total invariant mass squared of all clusters entering the jet algorithm. The 
definition of mf^ is taken to be 2EiEj (1 — cos 6ij). Here Ei and Ej are the energies of the 
clusters i and j, and 9ij is the angle between them. 

In its conventional form, the jet algorithm combines the pair of clusters i, j with the min- 
imum m'^j/W'^ to be a 'proto' jet by adding the four-momenta pi and pj. This prescription is 
repeated iteratively for the remaining clusters and 'proto' jets until all possible combinations 
i, j lead to rn^j/W'^ > ycuu the jet resolution parameter. In the present analysis, however, we 
use the jet algorithm to recombine the accepted clusters iteratively up to the point where ex- 
actly (2-1-1) jets remain. The smallest scaled jet mass given by any combination of the (2-1-1) 
jets is defined to be the observable ?/2- The y2 distribution, l/N^is dn/dy2, where Ncis is the 
number of deep-inelastic scattering events passing the selection of section 3, corresponds to the 
differential (2-1-1) jet event rate. 

The same definition of the (2-1-1) jets and of the variable y2 is used for the analysis of the data 
and of the Monte Carlo events after detector simulation. In events simulated at the hadron or 
parton level and in the NLO calculations, the jet algorithm is applied to hadron or parton four- 
momenta, respectively. The polar angle cut of 7° which is applied for clusters is also applied 
for hadrons, 6had, and partons, 9par- We take all components of the 'missing momentum' due to 
this cut into account and do not neglect the mass. 

With these definitions we observe that the smallest mass rriij of all possible combinations 
i,j of the (2-1-1) jets is most likely to be obtained by the combination of the two non-remnant 
jets. The fraction of events in which the minimum mass is formed by inclusion of the remnant 
jet is of the order of 15% for both data and NLO calculations. 

In addition to the definition given above, we also measure the differential (2-1-1) jet rate and 
determine «s using the E, Eq and P variants of the JADE algorithm [ [TtI ] without performing a 
full analysis of systematic errors. For these three algorithms mf^ is defined as (p, + Pj^. For 
the E algorithm, the combined four-momentum is simply the sum of the four-momenta pj 
For the Eq algorithm, the combined energy is defined as Ei + Ej and the combined momentum 
is ^^f'l (Pi + Pj)- For algorithm, the combined momentum is pi, + pj and the combined 
energy is \pi + Pj\- The definition of the recombination scheme for the latter two algorithms 
implies that the reconstructed jets are massless. This is not the case for the JADE algorithm and 
its E variant which conserve energy and momentum exactly in the recombination procedure. 



5 Description of the data by QCD models 

Before correcting the (2-1-1) jet event rate for detector and hadronization effects, as described 
in the next section, we study the description of the data by the QCD models LEPTO 6.5 and 
ARIADNE 4.08. LEPTO is based on the exact first order matrix elements followed by higher 
order radiation approximated by leading logarithm parton showers. In contrast, ARIADNE 
models the QCD cascade by emitting gluons from a chain of radiating colour dipoles. In QCD- 
Compton events the dipole is formed between the struck quark and the proton remnant, and 
the first gluon emission reproduces the first order matrix elements. In boson-gluon-fusion 
events, the quark and the antiquark are generated according to the first order matrix elements. 
Two dipoles are formed between each quark and the proton remnant and continue to radiate 



independently. Both LEPTO and ARIADNE use the Lund string hadronization model Qlq]. 
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We used the parameters of LEPTO and ARIADNE tuned to reproduce published HERA data 
[[T9|], in combination with the parton density functions of MRSH [pO|]. The generated events 
were passed through a full simulation of the HI detector. For each model an event sample was 
generated that was ~ 6 times larger than that of the experimental data. The same event and 
cluster cuts are applied to the simulated events as to the data. 

In Figure |l] we show the distributions of four representative jet variables: the differential 
(2+1) jet event rate y2, the variables Zp and Xp, and the polar angle of the most forward jet. The 
variable y2 was defined above. The definitions of Zp and 

Zj, = min Ei (1 — cos 6i) / > Ei(l — cos 9i) and x-n = — ^ ^ , 

where Ei and 6i are the energies and polar angles of the two non-remnant jets remaining af- 
ter the clustering of the jet algorithm, and is the corresponding invariant jet mass calcu- 
lated without neglecting the jets' masses. The variables Xp and Zp measure the approach to 
the (2-1-1) (1-1-1) singularities corresponding to the two non-remnant jets becoming one jet 
(Xp — ^ 1) or as one jet is absorbed into the remnant jet {Zp 0). 

In order to increase the fraction of events with a clear (2-1-1) jet structure, thus enhancing the 
sensitivity to ag, we define a subsample of events with y2 > 0.01. To decrease the sensitivity 
to the modeling of initial-state multi-parton emissions and to avoid forward jets which are too 
close to the proton remnant, we require that the non-remnant jets satisfy 10° < 9 jet < 145°. 
The requirement 9jet > 10° is found, in particular, to improve the description of the data by 
LEPTO. After these cuts, the (2-1-1) jet event sample consists of 2 235 events. 

In Figure |T](a) the uncorrected y2 data distribution is compared with the predictions of 
LEPTO and ARIADNE. The distribution is normalized to the total number of DIS events Nois 
selected in section 3. Both models give an acceptable description of the data. At large values 
of y2, the distribution from ARIADNE tends to be above that of the data while LEPTO is sys- 
tematically low. In Figure |l](b) and (c) the Zp and Xp distribution are shown for uncorrected data 
and the models mentioned above. ARIADNE roughly describes the measured Zp distribution 
with the exception of the first bin, while LEPTO and data disagree in particular in the lowest 
two Zp bins. The poorest description of the data is observed for the Xp distribution. ARIADNE 
approximately reproduces the data in the central part of the distribution. It overestimates and 
underestimates the data in the very low and very high Xp region, respectively. LEPTO shows the 
opposite trend. Note that the drop of the Zp distribution at the lowest Zp bin and the decrease of 
the Xp distribution at large values of Xp are consequences of the cut y2 > 0.01. The distribution 
of the polar angle of the most forward non-remnant jet is shown in Figure |T](d). It is sharply 
peaked at small angles and is well described by both models. 

We have studied the accuracy with which the data is described by ARIADNE and LEPTO 
for a wide range of selection criteria in addition to those discussed above. Overall ARIADNE 
gives the better description of the data. We conclude that the qualitative description of the data 
is acceptable and that a one-dimensional correction of the y2 distribution is possible although 
an improved model description of the data is clearly desirable. In the following analysis, we 
correct the measured y2 distribution with ARIADNE and use LEPTO as a consistency check. 

6 Correction of the data 

We correct the measured 2/2 distribution by the method of regularized unfolding described in 
[0. First, we unfold the y2 distribution for detector effects only, in order to make direct 
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comparisons with QCD model predictions possible. For each simulated ARIADNE event, the 
value of y2 is determined by clustering hadrons and simulated calorimeter clusters, respectively. 
Then, the y2 distribution calculated from hadrons is reweighted such that the y2 distribution 
from simulated clusters best fits the data. The weights are found by means of a log-likelihood 
method where strongly oscillating solutions are suppressed. As result, we obtain four bins of a 
reweighted y2 distribution - corresponding to unfolded data. The unfolded distribution is given 
in Table [T]. The quoted systematic error consists of two contributions added in quadrature: the 
influence of the uncertainty of the absolute hadronic energy scale of the liquid argon calorime- 
ter, and the full difference to the y2 distribution unfolded with LEPTO instead of ARIADNE. 
The unfolded y2 distribution is shown in Figure ^ together with the predictions of LEPTO and 
ARIADNE. The statistical error is of the order of 5% but the systematic errors can be larger. 
Both models roughly reproduce the data. The prediction of ARIADNE is high at large y2, 
while LEPTO falls too low. These observations are consistent with our conclusions from the 
comparison of the uncorrected data and the predictions of LEPTO and ARIADNE in Figure |l]. 

Next, we unfold both detector and hadronization effects in a one-step procedure in order to 
compare the data to NLO predictions. The unfolded distribution is also listed in Table |T] and is 
discussed in the next section. The size of the combined hadronization and detector migration is 
illustrated in Figure ^ where the reconstructed y'^^^ after hadronization and detector simulation 
is compared with the yf"'^ found by clustering the partons before hadronization. The bins shown 
correspond to those selected for the determination of as- For both LEPTO and ARIADNE a 
significant correlation is observed between y^^^ and yf""^. The y2 distribution is systematically 
shifted to smaller y2 values after hadronization and detector simulation, and the migrations are 
sizable. This is why a full unfolding procedure is used as opposed to a bin-by-bin correction 
factor method. We study the systematic uncertainty of the migrations in detail in section 8 by 
using alternative QCD models for the correction of the data and by varying model parameters. 
We also compared distributions of other jet variables like Xp, Zp and jet polar angles for partons 
and for reconstructed clusters after hadronization and detector simulation. All the jet variables 
show clear correlations between the different levels. 



7 NLO predictions and determination of a 5 
7.1 NLO QCD programs 

The NLO predictions are calculated with MEPJET, version 1 .4 [|Sp . MEPJET allows arbitrary 
jet definitions and the application of cuts in terms of parton four-momenta. Other programs 
were limited to a specific jet algorithm and made approximations in regions of phase space 
relevant for previous as analyses ^] that turned out to be imprecise [^. MEPJET uses a 



'phase space slicing' method p4| ] to deal with final-state infrared and coUinear divergences 
associated with real emissions of partons. If the invariant mass squared s of a pair of partons in a 
multi-parton state is smaller than a technical parameter Smin, soft and coUinear approximations 
are applied to perform the phase space integrations analytically. The infrared and coUinear 
divergences thus extracted cancel against those from the virtual corrections. If s exceeds Smin 
the integrations are done numerically without using explicit approximations. 

We run MEPJET with Smin set to the recommended value of 0.1 GeV^. The statistical 
precision of the predicted y2 distribution is ~ 1%. As a cross check, we changed Smin from 0.1 
to 0.05 and 0.01 GeV^ in MEPJET and observed no significant changes in the y2 distributions. 
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Note that our statistical precision at Smin = 0.01 GeV^ is then reduced to ~ 2% due to the larger 
fraction of (3+1) parton states treated numerically. 

More recently the program DISENT ^ became available which uses a different technique 
to treat divergences based on a 'subtraction' method p3] in combination with dipole factoriza- 



tion theorems [pq]. While we use MEPJET for this analysis, we have compared the predictions 



of MEPJET and DISENT version 0. 1 , which were run with the same value of A^IL. and the same 

' MS 

parton densities, for all crucial distributions of this analysis and find general agreement at the 
level of a few percent. Looking to the y2 distributions in detail, however, we see a significant 
discrepancy which is of little relevance for this analysis and is translated into an error in in 
section 8. To leading order we find the predictions of MEPJET and DISENT to be consistent 
within a fraction of a percent. 



7.2 Comparison of QCD model and NLO predictions 

Before extracting a value of from a comparison of corrected data and NLO calculations, a 
region of jet phase space must be identified in which NLO predictions can provide a fair de- 
scription of jet related observables. We verify the extent to which this is the case for the above 
selections by comparing NLO jet distributions with the parton jet distributions of ARIADNE 
and LEPTO. The use of QCD model predictions rather than corrected data distributions gives 
reduced statistical error. In addition, the comparison of ARIADNE and LEPTO provides inter- 
esting information on possible ambiguities in the definition of the parton level to which the data 
are corrected. 

In Figure ^a) the ?/2 distributions for ARIADNE and LEPTO are shown together with NLO 
calculations for different values of A^i^. In order to avoid a dependence of the following study 

on the value of a, we chose the extreme values of A^^ = 100 MeV and 600 MeV corre- 

^ MS 

sponding to q;s(M|) = 0.097 and 0.132, respectively. (Note that A^IL. serves only as a technical 
steering parameter for MEPJET.) The number of flavours used in the calculation is set at five. 
As with ARIADNE and LEPTO, the MRSH parton density functions are used in MEPJET. The 
same cuts on the hadronic final state, y2 > 0.01 and 10° < Ojet < 145°, that were applied for 
Figure 1 are used here. Note that the mean number of partons per event with 9par > 7° is 9.7 
for ARIADNE and 10.7 for LEPTO, whereas in MEPJET at most three partons and the proton 
remnant are produced. 

We find that the distributions derived from ARIADNE and LEPTO are in qualitative agree- 
ment, and that their shapes are similar to those of the NLO distributions. However at larger val- 
ues of y2, ARIADNE approaches the MEPJET prediction for A^L = 600 MeV, while LEPTO 

comes closer to the A^^ =100 MeV. This trend corresponds to that observed from the com- 
parison of data and ARIADNE and LEPTO after detector simulation as shown in Figure 1(a). 

The predictions of MEPJET and the distribution from ARIADNE for the Zp variable, shown 
in Figure ^(b), are also in fair agreement. LEPTO falls below ARIADNE at small Zp which 
is also seen in Figure [l|(b). A pronounced difference between ARIADNE and NLO is seen in 
the Xp distributions shown in ^c). This effect is not sensitive to changes of the phase space 
selection criteria and is further discussed in section 8. The corresponding prediction of LEPTO 
agrees well with that of the NLO calculations. The distributions of the forward jet's polar angle 
from ARIADNE and LEPTO which are shown in Figure ^d) are well described in shape by 
QCD in NLO. 
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We conclude from this comparison that within the phase space region selected by the cuts 
listed above, NLO calculations are expected to provide an adequate description of jet produc- 
tion in the data. This statement remains qualitative at this stage since we do not yet make an 
assumption on the value of as to be used in the NLO calculations and since we observe system- 
atic differences between ARIADNE and LEPTO. 



7.3 Fit of a > 



The y2 data distribution corrected for detector and hadronization effects is compared with MEP- 
JET in Figure ^. For the first and last bin in particular the systematic error (see section 8) is 
large compared with the statistical error and is dominated by the model dependence. The NLO 
predictions of MEPJET for different values of A^i^ are also shown. 

In NLO the differential jet rate is given by the expansion I/ctdis da2+i/dy2 = ^(^2) «s + 



MS 



100 and 



B{y2) a^. From the y2 distributions in NLO, obtained by running MEPJET for A. 
600 MeV, we obtain the coefficients A and B for the four bins in y2 evaluating at the scale 
= < Q"^ >, where < Q'^ >~ 620 GeV^ is the mean Q'^ of our (2-1-1) jet event sample. The 
mean of the entire selected DIS event sample is 545 GeV^. 

In order to relate A^ to A^ and thus to at a given scale u^, we use the following 



formulae [ETl ^ 



MS 



^^MS ^^MS y^^MS'"''^^ 



In(m2/A.(') 



MS 



-963/13225 



with nib, the mass of the bottom quark, set to 5 GeV, and the two-loop expansion 



An 



1 - 



2(3i lnln(/iVA 
^0 ln(/i7A 



MS 



MS 



19 



51 — ^ nj, and Uf the number of quarks of mass less than /i. 



with (3q = 11 — ^ Uf and /3i 

namely rij = 5 in our case. The same formulae are used in the MEPJET program. Given A and 
B and the relation of A^ to a^, the NLO j/2 distribution can conveniently be calculated for any 



value of A 



MS' 



We perform a minimum fit (x^/d.o.f. = 6.9/3) of A^^ taking into account the statis- 
tical correlations between the bins of the unfolded data distribution. As the result we obtain 
AgL = 320 ± 33 MeV corresponding to a, (Af|) = 0.118 ± 0.002 (stat.). Note that the 
choice of fx^ = < > for the calculation of the coefficients A and B is to some extent 
arbitrary. It influences the value of A and B but not the value of the fitted A^^. The NLO 



MS 



prediction corresponding to the fitted value of Ai^L is shown as the full line in Figure and a 
good description of the data is observed. 
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8 Determination of systematic errors 



We study various effects that might influence the result by varying the hadronic energy scale 
of the liquid argon calorimeter, changing the experimental cuts, and by using different Monte 
Carlo models for the data correction. We also use alternative parton density functions, measure 
jet rates with different variants of the modified JADE jet algorithm, and choose different renor- 
malization and factorization scales. The various fitted values of corresponding to different 
classes of uncertainties are shown in Figure ^ All values of given in the following refer to 
as at the scale /x^ = M|. 

Energy calibration 

The hadronic energy scale of the liquid argon calorimeter is varied by ±4% which leads to a 
systematic shift in ?/2- Note that there is no fully compensating effect in the ratio ?/2 = mfj/W'^ 
due to the definition of W which includes the 'missing momentum' vector. The resulting un- 
certainty in «s is ±0.003. 

The variation of the electromagnetic energy scale of ±3% leads to a negligible change in ag. 
Polar angle cuts Odus/Opar 

We vary the minimum value of the cluster acceptance cut Odus and in parallel the corresponding 
cut for partons 9par within a range of 5° — 15°. The variation of the Odus/dpar cut checks the 
quality of the detector simulation mostly but also the description of the data in the forward 
detector region where the models are less well tested. It also shows the stability of the proton 
remnant separation by the jet algorithm. 

The additional values fitted in this range of cluster cuts are slightly lower than the main 
value, the smallest one differing by —0.002. We see no indication for a systematic trend in ag as 
a function of the cut value. Without any cluster or parton cut the qualitative agreement between 
NLO and ARIADNE/LEPTO parton distributions deteriorates and stricter phase space cuts are 
needed. As an example we omit the cluster or parton cuts as well as the forward jets polar angle 
cut of 10° but apply the additional event cut Zp > 0.15. This reduces our (2+1) jet event sample 
by roughly a factor of 2. We obtain an value of 0.119 ± 0.003 (stat.) which is consistent 
with our main result. 

Event selection cuts 

In addition to the OduslQpar cut variation, we study the variation or introduction of various event 
cuts. As before, all cuts are applied in parallel to quantities calculated from measured clusters, 
from simulated clusters and partons, and from the partons of the NLO program. We require 
?/2 > 0.02 instead of 2/2 > 0.01. We change the polar angle jet acceptance cut to 6 jet > 8°, 
12° or 14°. We require Zp to be larger than 0.05, 0.1 or 0.15. We unfold the differential jet rate 
for different Q"^ ranges and vary the minimum cut from = 200 GeV^ to = 100 and 
250 GeV^. Most of these variations correspond to significant changes in the number of events 
considered. However we find a variation of ag of +0.002 and —0.003 at most. The largest 
change of —0.003 is found for the cut ?/2 > 0.02. Note that the sizes of the observed changes in 
as are close to those of our statistical error, and that no indication of any systematic trend as a 
function of a cut variation is visible. Thus we regard the analysis as stable with respect to the 
phase space selection. 
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Model dependence 

We test the model dependence of the result by repeating the analysis using LEPTO for the 
correction of the data. When using cuts identical to those given before, a value of = 0.1 16 
is obtained. This result is reasonably stable with respect to the variation of 9cius, Ojeu Zp and y2 
cuts, although the observed changes of the determined as values are larger than for the analysis 
based on ARIADNE. 

Motivated by both the poor agreement of the shape of the Xp distributions between ARI- 
ADNE and NLO in Figure ^(c) and the relatively large differences between the Xp distribution 
of data and ARIADNE in Figure [T](c), we reweighted ARIADNE events such that the measured 
Xp distribution is reproduced. Effectively, this can be seen as a correction to the parton evolution 
mechanism of ARIADNE. This procedure leads to negligible change in the corrected y-2 data 
distribution but we find better agreement between ARIADNE and NLO in Figure ^. Reweight- 
ing ARIADNE in Xp also gives a good description of the dependence of the rate of (2+1) jet 
events, i?2+i(<5^) = N2+i{Q'^)/NDisiQ'^), where ARIADNE (unweighted) was shown to be 
inferior to LEPTO 



Possibly large hadronization corrections could fake the radiation of hard partons described 
by perturbative QCD, and could cause systematic biases in the correction of the data. The un- 
certainty of the hadronization corrections is not directly tested by the comparison of LEPTO 
and ARIADNE since both models use the Lund string hadronization. We thus vary the param- 



eters a and b of the Lund fragmentation function [ [181 ] and the parameter ag, which determines 



the mean pt of a produced hadron, from their default values a = 0.3, b = 0.58 GeV ^ and cXg 
= 0.36 GeV, to a = 0.1 and 1.0, 6 = 0.44 and 0.70 GeV-^ and to = 0.25 and 0.45 GeV. We 
derive hadronization correction factors for the y2 distributions obtained from the events simu- 
lated with these sets of parameters. The differences between calculated correction factors do 
not exceed a few percent, and the corresponding variations in ag which we find are at most 
+0.002 and —0.002. In the same manner we vary the parameter Qq of LEPTO which cuts off 
the evolution of the final state parton shower. Setting Qq to 4 instead of 1 GeV we observe a 
change of +0.004 in as- Setting the corresponding parameter for the initial state parton shower 
from 1.5 to 4 GeV, we observe a change of only 0.001 in as- 



We repeat the analysis with the QCD model HERWIG, version 5.8 [Q. HERWIG com- 
bines a model for coherent parton shower radiation and an additional first order matrix element 
correction. Hadronization follows the cluster fragmentation model [Pl]]. The description of the 
data with HERWIG is satisfactory for our purpose although HERWIG does not describe the j/2 
distribution at very small y2 and predicts the fraction of (2-1-1) jet events to be ~ 10% lower 
than that of the data. Unfolding the data with HERWIG leads to a change in the fitted as value 
of —0.006. As result of the described variation of the models and of the model parameters we 
assess the total model dependence of our measurement to be +0.004 and —0.006. The model 
dependence represents the main source of experimental uncertainty. 



Parton density functions 

The fit to the experimental y2 distribution is repeated for several choices of parton density func- 
tions in MEPJET including GRV HO (92), CTEQ2pM and CTEQ4M. We find a maximum 
variation of +0.005 and of —0.001. This dependence is mostly due to the uncertainties in the 
gluon density function. Gluon-initiated processes account for ~ 50% of the (2-1-1) jet events in 
our sample. 

Since we run MEPJET for values of different from those assumed during the global 
fits to deep-inelastic scattering data in which the parton density functions were determined, we 
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study the effect of this inconsistency. This is done using the MRS A', MRSR and CTEQ4A 
series of parton density functions which each combine parton density functions determined on 
the basis of the same experimental data and the same fit procedure but with set to different 
values. From the observed change in the fitted value of «s we estimate the effects of this incon- 
sistency to be smaller than ±0.002. In total, we assign an uncertainty of +0.005 and —0.002 
due to the knowledge of the parton density functions. 

Different jet algorithms 

In addition to the JADE algorithm, we determine ag with three related cluster algorithms, 
namely the E, Eq and the P algorithms. The unfolded differential jet rate distributions are 
given in Table |^. 

Comparing the measured y2 distributions for the JADE, E and P algorithms, we observe 
small but statistically significant differences. Similar differences are observed for the corre- 
sponding NLO predictions, which are given in Table |^, such that the fitted values of ag do not 
differ much. We obtain a, (M|) = 0.119 ± 0.002 and a,(M|) = 0.117 ± 0.002 for the E 
and P algorithms, respectively. 

The measured y2 distributions for the Eq algorithm is closest to that for the JADE algo- 
rithm. To next-to-leading order these algorithms are identical. This is a consequence of the 
jet finders' definitions and of the fact that in this analysis no cuts on the jets' transverse (or 
longitudinal) momenta are made. The (small) difference in the measured y2 distributions from 
these two algorithms may be interpreted as an expression of higher order recombination effects 
which cannot be accounted for in 0(Q;f) calculations. The value of determined with the Eq 
algorithm is as(M|) = 0.120 ± 0.002. The observed differences between the results of the 
different algorithms are small and are not treated as an additional error. 

Renormalization and factorization scale 

In NLO the y2 distribution depends on the choice of the renormalization and factorization scales 

and /ij. We estimate the renormalization scale dependence by varying fi^ from to 1/4 
and 4 in MEPJET and by repeating the as fit. The corresponding uncertainty in as is +0.007 
and —0.005. In addition, we use the scalar sum of the transverse momenta of the jets in the 
hadronic centre-of-mass frame as a renormalization scale. This corresponds to a considerable 
difference in the magnitude of the renormalization scale given that is typically about a factor 
of 20 larger than the square of a jet's transverse momentum in the hadronic centre-of-mass 
frame. The corresponding change is close to that observed for /i^ = 1/4 Q^. 

In deep-inelastic scattering, coUinear initial-state divergences are absorbed into redefined 
parton densities introducing the dependence on a factorization scale. We estimate the factoriza- 
tion scale dependence by varying n'j from to 1/4 and 4 as we do for fi^, and we find 
the factorization scale dependence of as to be ±0.001. 

DISENT NLO predictions 

As an alternative to the determination of as based on MEPJET calculations we use DISENT 
predictions for the fit. The resulting change in as is —0.003. 

Discussion of higher order effects 

This and similar analyses of e+e^ annihilation data rely on the assumption that the distri- 
butions of observables obtained from NLO calculations and from the partons in the parton 
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shower/colour-dipole models are comparable. There are, however, ambiguities in the definition 
of the parton level of the QCD models used to correct the data which can lead to uncertain- 
ties in the determined value of a^. One may argue that the data should not be corrected to the 
level of the jets reconstructed from the final partons before hadronization but to that of the jets 



reconstructed from the partons at an earlier stage of the parton shower [ |33| ] 



We study the dependence on the correction levels in two different ways. First, we in- 
vestigate systematic changes of y2 due to subsequent recombinations of partons during the 
dipole/parton shower radiation. We 'pre'cluster the final partons of ARIADNE or LEPTO using 
mf,j = 2 min {Ef, Ej) (1 — cos %) as used in the Durham jet algorithm [ ]T7| ] and then continue 
clustering using the JADE definition. The 'pre' clustering is stopped when the scaled invariant 
masses of all pairs of parton jets satisfy m'^j/W'^ > 0.00005, where rriij is calculated accord- 
ing to the JADE definition. This corresponds to the point where on average (4-1-1) parton jets 
remain. We apply the JADE algorithm to these parton jets, continue clustering up to (2-1-1) jets 
and calculate y2. Comparing the result using this procedure with the y2 distribution obtained 
using the JADE algorithm throughout, we see differences of a few percent for either ARIADNE 
or LEPTO. This translates into a similar difference of a few percent in the result for a^. Larger 
differences are observed, however, when we extend the 'pre'clustering further. 

In the same spirit, we also change the value of the parameter Qq of LEPTO, which cuts off 
the final-state parton shower, as a means of looking at an early stage of the parton shower by 
forcing a change in the average number of partons produced. In contrast to the above variation 
of Qo in the context of the model dependence, here we are interested in the y2 spectrum of a 
(variable) parton level keeping the hadron level fixed. We compare the differences in the y2 
distributions obtained. Again changes of a few percent are observed for e.g. Qq = 3 GeV 
instead of Qo = 1 GeV, but differences increase with larger values of Qq. 

Note that the differences in the aforementioned definitions of the parton level may partially 
be due to unknown higher order corrections missing in NLO. The interpretation of these stud- 
ies becomes difficult and we have not included such estimates in our combined error. Such 
effects might, in principle, be detected by comparing the shape of the distributions from NLO 
calculations and from parton shower/colour-dipole models of Figure 0. From another view- 
point differences in the measured as values using different jet recombination procedures, or a 
large dependence on the renormalization scale could be symptoms of the same difficulty. The 
dependence on the renormalization scale is the dominant uncertainty which is included in the 
systematic error. This situation with higher order effects and/or the ambiguity of the parton 
level is similar to that in e^e annihilation where such uncertainties turned out to be an 
important limitation. 



Combined systematic error 

We define the combined systematic error in the following way: assuming the errors of the 
different classes to be largely independent of each other, the positive and negative systematic 
errors of the first 5 classes are each added in quadrature. Thus we determine the systematic 
error of this analysis to be +0.007 and —0.008. An additional theoretical error of +0.007 and 
—0.006 is obtained correspondingly, considering the measured difference between the JADE 
and Eq algorithm, the renormalization scale uncertainty and the uncertainty due to observed 
difference of MEPJET and DISENT. 
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9 Summary 



We have presented a measurement of jet related distributions in deep-inelastic scattering pro- 
cesses at HERA in the kinematic range 200 < < 10000 GeV^. The jets are found with the 
modified JADE jet algorithm. The measured jet distributions are compared with QCD model 
expectations and for most distributions we find acceptable agreement between the data and the 
models ARIADNE 4.08 and LEPTO 6.5. Acceptable agreement is also observed for the differ- 
ential (2-1-1) jet event rate after correcting for detector effects. 

The differential jet rate, corrected for both detector and hadronization effects, is compared 
with NLO QCD calculations in a region of jet phase space where the effect of higher-order 
parton emissions not considered in NLO is estimated to be small. A fit of the NLO predictions 
as a function of the strong coupling constant is performed which results in 



a,(M|) = 0.118 ± 0.002 (stat.) i^yst.) t'^oml (theory). 



A good description of the corrected differential jet rate by the next-to-leading order pre- 
diction is observed for the fitted value of a^. The resulting as is compatible with previous as 
determinations based on the same observable in e+e^ annihilation [|5|1 and with the world aver- 



age value of as [ p8[ | which provides a direct consistency check of perturbative QCD. The same 
conclusions are reached considering the results obtained with the E, Eq and P variants of the 
modified JADE algorithm. 

The most important uncertainties of the as value determined are caused by the as yet limited 
precision of the data description by current QCD Monte Carlo models, by ambiguities in the 
definition of the parton level to which the data are corrected, and by the large renormalization 
scale dependence. 
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y2 


1/Ndis dn/dyl^'"^ 


^/Ndis dn/dyf" 


d(T2+l/dy2 


0.010 - 0.020 
0.020 - 0.035 
0.035 - 0.055 
0.055-0.100 


13.43 ±0.491^:^^ 

4.67 ± 0.20 IgJI 
1.51 ± 0.08 

0.39 ± 0.02 


15.06 ± 0.87 1^:^^ 
5.14 ± 0.32 
1.71 ±0.12;°:^^ 

0.47 ±0.03 


13.70 ±0.15 
5.01 ± 0.06 
1.79 ±0.03 
0.53 ±0.01 



Table 1: y2 distribution determined with the JADE algorithm corrected for detector effects 
(I/Ndjs dn/di/2"''), corrected for both detector and hadronization effects (I/Ndjs dn/d|/f"'), 
and the NLO prediction obtained from MEPJET for as{Mz) — 0.118 in combination with the 
parton density functions MRSH {I/gdis d(T2+i/dy2)- All distributions are determined in the 
kinematic region defined in section 3 and the cut 10° < 6 jet < 145° is applied for hadron 
and parton jets, respectively. The first error is statistical, the second systematic. For the NLO 
calculation only the statistical error is given. 



par 


E 


£"0 


P 


0.010 - 0.020 


14.97 ±0.79 


15.14 ±0.84 


13.00 ±0.76 


0.020 - 0.035 


6.59 ±0.31 


5.54 ±0.32 


4.66 ±0.27 


0.035 - 0.055 


2.16 ±0.13 


1.76 ±0.12 


1.51 ±0.10 


0.055-0.100 


0.63 ±0.05 


0.49 ±0.03 


0.42 ±0.03 



Table 2: y2 distribution corrected for both hadronization and detector effects for the E, Eq and 
P algorithms. The error is statistical only. 



y^^O 


E 




P 


0.010-0.020 


15.59 ±0.15 


14.27 ±0.16 


12.39 ±0.15 


0.020 - 0.035 


5.67 ±0.06 


5.13 ±0.06 


4.42 ± 0.06 


0.035 - 0.055 


2.10 ±0.03 


1.85 ±0.03 


1.53 ±0.03 


0.055-0.100 


0.66 ±0.01 


0.53 ±0.01 


0.46 ±0.01 



Table 3: NLO y2 distribution obtained from MEPJET for the E, Eq and P algorithms for the 
fitted values of ^^(Ml) — 0.119, 0.120 and 0.117, respectively, in combination with the parton 
density functions MRSH. 
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Figure 1: Distributions of (a) y2, (b) Zp, (c) Xp and (d) the forward jet's polar angle Of^djet for 
uncorrected data compared with the prediction of ARIADNE 4.08 and LEPTO 6.5 including 
full detector simulation. For all distributions the (2+1) jet event cuts y2 > 0.01 and 10° < 
9jet < 145° are applied with the exception of (a) where the cut 1/2 > 0.01 is omitted. The 
distributions are normalized to the number of deep-inelastic events N^js passing the kinematic 
cuts. The errors are statistical only. 
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Figure 2: Distribution of the observable y2 corrected for detector effects compared with the 
prediction of the models ARIADNE 4.08 and LEPTO 6.5. As in Figure 1, the cut 10° < Oj^t 
< 145° is applied. The error bars correspond to the statistical and systematic errors added in 
quadrature. The inner error bars give the statistical error only. 
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Figure 3: (a) The distribution of ^2^"^ as reconstructed from the calorimeter clusters after 

hadronization and detector simulation from the events with 0.01 < yf"'^ < 0.02 as predicted by 
ARIADNE 4.08 (white circles) and LEPTO 6.5 (full circles). The distribution is normalized to 
the number of events with y^'^'^' in the range 0.01 < yf"*^ < 0-02. Figures (b), (c) and (d) show 
the same for different ranges of indicated by the legend and by the dashed vertical lines. 
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Figure 4: NLO predictions based on MEPJET for the distribution of (a) 1/2, (b) z,^, (c) x.^ and of 
(d) Ofwdjet compared with parton jet distributions of ARIADNE 4.08 and LEPTO 6.5, respec- 
tively. The full line corresponds to MEPJET for the extreme value of = 600 MeV and the 

dashed line to = 100 MeV. The cuts y2 > 0.01 and 10° < Ojet < 145° were applied for 
MEPJET, ARIADNE 4.08 and LEPTO 6.5 each. 
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Figure 5: Distribution of the differential jet rate y2 corrected for detector and hadronization 
effects compared with the NLO prediction of MEPJET for = 600 MeV (dotted Une) and 

= 100 MeV (dashed Une). The full line shows the NLO prediction for the fitted value 

of as which corresponds to A^ = 320 MeV. The error bars on the corrected data distribution 
correspond to the statistical and systematic errors added in quadrature. The inner error bars give 
the statistical error only. 
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Figure 6: List of systematic uncertainties on the fitted value of ag. The horizontal lines separate 
different classes of uncertainties. The vertical line indicates the central value of as{Mz) result- 
ing from the fit. The black points give the values of as that are obtained when each source of 
sytematic error is varied as described in the text. The uncertainties of the three classes below 
the dashed horizontal line are combined to give the theoretical error. 



25 



